Random circuits by measurements on weighted graph states 



OO 
O 

o 

(N 
G 

OO 



43 ■ 

Oh! 
-f— > 

cr 



> 

OO 

in 
o 

m 

o 

oo 
o 



X 



A. Douglas K. Plato/'H Oscar C. Dahlsten, 2 'Q and Martin B. Plcnio 1 ^ 

1 Institute for Mathematical Sciences, Imperial College London, Prince's Gate, London SW7 2PG, UK 
QOLS, Blackett Laboratory, Imperial College London, Prince Consort Road, London SW7 2BW, UK 
^Institute for Theoretical Physics, ETH Zurich, 8093 Zurich, Switzerland 
(Dated: June 18, 2008) 

Random quantum circuits take an input quantum state and randomize it. This is a task with 
a growing number of identified uses in quantum information processing. We suggest a scheme 
to implement random circuits in a weighted graph state. The input state is entangled with the 
weighted graph state and a random circuit is implemented when the experimenter performs local 
measurements in one fixed basis only. The scheme uses no classical random numbers and is a new 
and natural application of weighted graph states. 

PACS numbers: 03.67.Bg, 03.67.Mn, 05.40.-a 



I. INTRODUCTION 

Randomising inputs has a wide usage in information 
processing. In the case of quantum information pro- 
cessing, the input may correspond to a pure state and 
the randomisation a unitary evolution picked at random. 
The random unitary cannot be generated directly, but in- 
stead needs to be implemented as a sequence of random 
elementary gates, i.e. a random circuit. 

The number of applications for random circuits in 
quantum information process is growing. They may be 
applied to hide information about input states to 
sample state spaces 0, and to model random processes 
such as thermalisati on p , d, [f| @, 0, H, 0| and black hole 
information leakage [id ]. A more surprising application 
is in the superdense coding of quantum states [ll( • Gen- 
erating uniformly random operators over U(2 N ) is, how- 
ever, exponentially hard. Instead it is often sufficient 
to produce circuits which are identical in only a few of 
the relevant statistical properties. Such pseudo-random 
circuits have already been efficiently implemented in an 
iVMi?-setup [13 • 

In the present work we propose a new experimental 
implementation for such random circuits. One prepares, 
in advance, a highly entangled resource state called a 
weighted graph state. The input state is entangled with 
a number of qubits in the resource state and the exper- 
imenter then performs local measurements in one fixed 
basis only. The randomness of the measurement out- 
comes effectively pick the circuit, and the output state is 
then carried by a subset of qubits of the total system. 

We identify several advantages of our scheme. No clas- 
sical random numbers are used in the process, so the ran- 
domness is entirely quantum. Furthermore, the measure- 
ment based approach has the advantage of being com- 
paratively scalable. The fact that the experimenter only 
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needs to measure in one basis throughout is a simplifi- 
cation and is to our knowledge at present unique to our 
scheme. The scheme is furthermore a new application of 
weighted graph states which can help to motivate such 
experiments. 

Recent work has independently suggested the measure- 
ment based approach for generating random circuits [l3| . 
They were able to improve on a previously known gate 
model, by first translating it into the conventional clus- 
ter state formalism. Here, however, we find that using 
weighted graph states is arguably more natural and sim- 
ple: no classical pseudo-random numbers are required, 
and measurements are only performed in one basis. 

We proceed as follows. Firstly background material on 
random circuits, measurement based quantum computa- 
tion and weighted graph states is presented. We then de- 
fine and analyse our proposed method. In the subsequent 
sections the scheme is evaluated, and finally we summa- 
rize possible experimental implementations and discuss 
open questions. 



II. BACKGROUND 
A. Random Circuits 

In this section we define random quantum circuits and 
then list properties we deem desirable, motivated by the 
consideration of possible applications. 

A sequence of quantum gates, i.e. unitaries, is called a 
circuit. Here a sequence of quantum gates picked at ran- 
dom is called a random circuit. Thus according to this 
definition, a simple example of a random circuit for a sin- 
gle qubit is the following: apply the Hadamard[36j] gate 
H, Pauli I or 1 each with probability 1/3. However, 
in practice, we are often more interested in random cir- 
cuits with particular properties, and so we include three 
additional conditions that such proposals should satisfy. 
These will serve as the basis for our later evaluation. 

(i) Unbiased sampling asymptotically 
For hiding an input state fully one requires the output 



states to have a flat distribution in state space, i.e. ac- 
cording to the unitarily invariant Haar measure P where 
P{\^)) = P(U\^)). It is a simple but important observa- 
tion that if the gates are picked from a non-universal set 
one would not access all states, and therefore not achieve 
the uniform distribution, even in the asymptotic limit. 
In fact it is essentially sufficient for a random circuit to 
consist of gates picked from a universal set of gates for it 
to induce the Haar measure asymptotically - see [l4| for 
discussion. 

(ii) Good sampling in polynomial time 
Unbiased sampling of the uniform distribution to within a 
fixed accuracy requires exp(N) elementary gates, where 
N is the number of qubits [h|. However for practical 
applications one will require the circuit to give a suf- 
ficiently good sampling in a modest, i.e. poly(N) time. 
What merits being called a sufficiently good sampling can 
depend on the application. For example, in supcrdense 
coding of quantum states 0, HH, [HI one requires that the 
entanglement is typically maximal [H, [l?], E, Oil • In 
fact, typically maximal entanglement may be achieved to 
some accuracy within poly(N) elementary gates picked 
at random. Such schemes exist [H, Q , results which were 
recently generalised in [2ll ], 

(Hi) Feasible experimental implementation 
For practical applications we furthermore require a feasi- 
ble experimental implementation, ideally one that is pos- 
sible with current technology. 

It should be noted that the unbiased sampling of prop- 
erty (i) is often included in definitions of random (and 
pseudo-random) circuits found elsewhere in the litera- 
ture. 



B. Measurement based quantum computation 

In this section we briefly introduce one-way, also called 
measurement based, quantum computing together with 
some relevant notation. Firstly we describe how to evolve 
a state by measurements followed by an introduction to 
the idea of resource states for measurement based quan- 
tum computation. For a more comprehensive introduc- 
tion we refer the reader to 
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A quantum computation is an evolution of an input 
state to an output state. The unitary transformation can 
in fact be realized through measurements. For example, 
say that one wants to perform the Hadamard gate, H, on 
a single qubit in state |\&) = a|0) We start by first 

preparing a second qubit in state |+) = -4= (|0) + |1)). 
Then we perform a Control-Z gate (CZ)[37| on the two 
qubits, yielding the state C Z (a\0) + (3\1})\+) = a|0)|+> + 
). Finally, we make a projective measurement in 
the |+)/|— ) basis on the first qubit. Now, if an outcome 
of +1 is obtained, the new state will be |+)(a|+) )). 
Thus the second qubit is now in the state H\^>), as de- 
sired. One notes that there is a random choice of the 
gate applied, depending on the two possible measurement 
outcomes. For controlled computations this randomness 



FIG. 1: A resource state for quantum computation on a single 
qubit together with the input state. Circles are qubits and 
the qubit (ini) carries the input state. The figure shows the 
state before the measurements have been performed. Solid 
lines indicate CZ gates. The resource qubits are in the |+) 
state before the CZ is applied. They are labelled m, ri2— rn, 
where I signifies the length of the circuit and the first and 
second index is the row and column in the graph respectively. 



is in general a problem, albeit a surmountable one [221 ]. 
However, for the purpose of producing random circuits 
this effect can, on the contrary, be an advantage. 

In order to perform a sequence of gates one can reit- 
erate this procedure, i.e. introduce a third qubit in state 
|+), apply CZ on qubits 2 and 3 and so forth. However 
the CZ operations commute with the relevant measure- 
ments in such a way that one can equivalcntly perform 
the CZ gates first, before doing any measurements. All 
the non-input qubits can be prepared in this manner and 
be thought of as a resource state, to be used together 
with an arbitrary input state. Figure [TJ describes a pos- 
sible resource state for the case of one qubit. This can 
naturally be generalised to more qubits, which will be 
discussed in the next section. 

It is often argued that the resource state measurement- 
based approach is easier to implement in experiment than 
the circuit model because one can do most of the entan- 
gling operations before the input state has been given. 
Since these gates are particularly difficult to do, it is de- 
sirable to be able to fail and redo them, without losing 
the input state. Only when this has been successfully 
achieved does one give the input state. This requires 
entangling the qubits carrying that state with those in 
the resource. After which, only local measurements on 
the system are necessary to achieve an arbitrary unitary 
transformation 12211. 



C. Weighted graph states 

There are many types of resource states one can use 
for measurement based quantum computation. In this 
section we define a type of resource state called a graph 
state [12]. Then we introduce a natural generalization 
thereof called weighted graph states, which is what we 
will use for our scheme. 

We begin by defining a graph, G = (V, E), finite 
non-empty point set V along with a collection E C V 
of unordered pairs of points in V. We say V is the col- 
lection of the points or vertices of G and E the collec- 
tion of edges. A graph is simple if it is undirected, and 
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has no loops or weights. We can describe a graph by 
the elements of a matrix T ab , where entries correspond 
to whether an edge connects vertices a and b. For un- 
weighted graphs, r a {, = 1 for all {a, b} <E E and zero 
otherwise. Such a matrix is called an adjacency matrix 
of the graph G, and so a simple graph can be described 
by a symmetric matrix, of ones and zeros, having zeros 
along the main diagonal. 

Graphs such as these can be used to describe a fam- 
ily of quantum states in the following manner. We first 
consider each vertex as labeling a qubit in the state 
|+). Between each pair of qubits connected by an edge 
in the associated graph we apply a unitary operation 
U ab (T) = e- ir <"*f( J -O®( J -0, i.e. a CZ gate. The 
resulting states described by this procedure are graph 
states, and includes the cluster state resources as spe- 
cial cases. Using this formulation, it is then straightfor- 
ward to generalise graph states to weighted graph states 
(WGS) , by simply relaxing the condition that edges carry 
no weights. Thus, the simplest definition of a weighted 
graph state is, 

\*wgs}= n u ab (r ab )\+)® N , (i) 

{a.b}eE 

where the product is taken over all edges {a, b} and the 
unitaries are now defined by U ab = e~ ira6 ' ! r( / ~ cr ")®( / - , 'z) ; 
where T ab are given by the components of the adjacency 
matrix of a weighted graph. 

Various applications and generalisations of weighted 
graph states themselves have been proposed 0, [25| . 
For example, one can include more general initial states, 
additional filtering operations and local unitary opera- 
tors. However, for the purposes of this paper, the above 
definition captures the relevant ideas needed to extend 
the cluster state formalism. 



III. MEASUREMENT BASED SCHEME 

We now describe our scheme for generating random 
circuits, using ideas from measurement based quantum 
computation. The key feature in this proposal is the 
adoption of weighted graph states as the resource states 
for measurement based quantum computing. This en- 
ables us to present a simple, fixed measurement algorithm 
which does not have the requirement of a separate clas- 
sical (or quantum) random number generator. It follows 
closely the simplest cluster state model, with only mi- 
nor variations in the non-local operations made between 
adjacent rows. 

First, we describe the proposed scheme. Then we char- 
acterise the associated evolution of the input state. This 
is important for showing that the process is unitary. 

The proposal is as follows. First, a planar rectangular 
array of N x I qubits is prepared, each initially in the 
state \ip} r - k = |+) with subscripts j and k labeling the 
row and column positions along the array, see FigfS] This 
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FIG. 2: The initial weighted graph state resource used in our 
random circuit generation scheme prior to any measurements, 
together with the column of qubits carrying the input state. 
Circles represent qubits and the qubits (im, ..Atin) carry the 
input state. Thin horizontal lines indicate CZ gates, while 
thicker vertical lines are <f> gates, defined in equation [2] The 
resource qubits are in the |+) state before the CZ and <j> gates 
are applied. Each lie on a vertex a — (j, k) in a graph G and 
take the labeling rjk. j is the row index which runs from 1 to 
N, the number of qubits in the input state, k is the column 
index which runs from 1 to I, the length of the circuit. 

configuration is represented by a graph G embedded in a 
2D lattice, and so a vertex, a, is described by two num- 
bers (J, k) 7 and on each vertex lies a qubit rjk- Between 
each pair of adjacent qubits on neighbouring columns we 
apply CZ gates, that is, U r]kTjk±l = diag(l, 1, 1, -1). 
Between each pair of adjacent qubits on neighbouring 
rows we apply the control phase gates, 

U(A jk r s±lk ) =diag(l,l,l,e-^™*), (2) 

where the parameters 4> ab = T ab j = <fi for all a, b. We will 
refer to this as the 0-gate. As these operators commute, 
they can be carried out in parallel. The number of rows 
needed is determined by the number of qubits on which 
the random unitary will act, while the depth or length / 
specifies the number of iterations to be performed. Note 
that it will not be necessary to prepare the entire resource 
state to depth I initially - it may be grown as the protocol 
proceeds. The circuit acts on an N qubit state \^)i n 
carried by qubits im, in% ■■■ in^. A CZ gate is applied 
between each irij and rji qubit. Then each in qubit is 
measured in the basis {|+), |— )}. The possible outcomes 
arc labeled by a bi-vector where each component 
Sj takes the value Sj = for a measurement yielding 

the eigenvalue 1 and Sj = 1 for the eigenvalue — 1. 

The measurements on the input column projects the 
first column of resource qubits into the state, 

|V>>®i=G«M(5«)|*> in (3) 
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where the operators acting horizontally between rows are 
given by 

M(S^) = Mi(s[ k} ) <g> M 2 (Si k) ) g> ■ ■ ■ <g> M n (S^), (4) 
with 

M j {sf ) )=HZ s ?\ (5) 
and the vertical operators by 



{( 3 ,k),(j',k)}eE 



(6) 



Where, again, the set of edges E is represented in Fig. 
[5J Measurements are then performed successively on 
columns 1 through l — l, leaving the output state on the 
hnal line of qubits. Immediately we can see that the evo- 
lution from column to column is unitary, and so the entire 
process can be described by a single unitary operation. 
In fact, as there is no feed- forward of the measurement 
outcomes one could in principle make all measurements 
simultaneously 

In Appendix [X] we consider more general weighted 
graph states and more general measurements. However, 
as the next section will show, the simpler version of the 
scheme appears to work very well. 



IV. EVALUATION OF SCHEME 
A. Unbiased sampling asymptotically 

We demonstrate numerically that the scheme gives un- 
biased sampling (Haar measure) in the asymptotic time 
limit. More precisely we show that the output state dis- 
tribution passes a necessary and stringent test, namely 
that the entanglement probability distribution approach 
that associated with the Haar distribution on pure states. 

The entanglement distribution associated with the 
Haar distribution has received considerable attention re- 
cently [IB, The average entropy of entan- 
glement TESa(Na, Nb) of a set of Na spins was studied 
already in the 70s and 80s [13, EU and the explicit so- 
lution ('Page's conjecture') was conjectured in [lj| and 
proven in [20( : 



MSa(N a ,N b ) 



1 
m~2 



2 N A + N 

E 

k=2 N B- 



•>N A 



1 



2-Wb+i 



(7) 



with the convention that Na < Nb and where Na + 
Nb = N, the total number of particles. 

This can be used to show that the average entangle- 
ment is very nearly maximal, meaning close to Na, for 
large quantum systems, i.e. N ^> 1. Hence one concludes 
that a randomly chosen state will be nearly maximally 
entangled with a large probability. Indeed, it was re- 
cently shown that the probability that a randomly chosen 




FIG. 3: (Color online) The average entanglement, E(S(pa)), 
(blue diamonds) of our randomized circuits compared with 
the average expected entanglement given by Page's conjec- 
ture (dashed line) for N — 8 and various Na- The solid red 
line shows the same quantities for stabilizer gates chosen at 
random. To calculate our averages we first perform 10 6 iter- 
ations on a randomly chosen initial state. Then we average 
the entanglement calculated over the next 10 6 iterations. In 
each case the fractional difference is less than 10 4 from that 
expected from true random circuits. 



state will have an entanglement Sa that deviates by more 
than 5 from the mean value W,Sa(Na, Nb) decreases ex- 
ponentially with S 2 [HI . 

Achieving this entanglement probability distribution is 
a strong condition to claim one has the true Haar distri- 
bution. Often one quantifies multipartite entanglement 
using "linearised" measures based on purity. For exam- 
ple, a popular measure used for pseudo-random circuits 
is the Meyer- Wallach entanglement [2(| , which is related 
to the average purity of each qubit. However, stabilizer 
states give exactly the correct purity one would expect 
from the Haar measure. Nevertheless they do not yield 
the correct entanglement distribution [27], [28| . The prob- 
ability distribution of entanglement, P(Sa) associated 
with stabilizer states sampled uniformly at random is 
given by Theorem I in [28] , 



nti(2 j + i) 



P(Sa) 



n 



N 

k=N-N A + 



i(2 fc + l) 



n ^ — i ( g ) 



2 2 i - 1 



where Na is the number of qubits belonging to Alice and 
N is the total number of qubits. Again Na < N ~ Na = 
Nb- The total state is bipartite and pure. 

Figures [3] and H] show that the output states do indeed 
follow the correct entanglement statistics. We also add 
the relevant stabilizer statistics for comparison, which 
highlight the fact that averaging over stabilizer states is 
not as good in this regard. 

The above considerations strongly indicate that the 
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states coming out of the scheme in the asymptotic time 
limit are indeed Haar distributed, but they do not con- 
stitute an analytical proof [Hj]. 



B. Good sampling in poly(N) time 

In the previous section we considered the sampling in 
the infinite time limit. In practise one will have a finite 
time available, and this time should scale as a polynomial 
function in N, the number of qubits carrying the input 
state. Under this restriction one cannot generate the 
Haar measure exactly, however one can hope to obtain a 
sufficiently good sampling of the state space nevertheless. 

To test this for our current scheme, we again use the 
entanglement of the output states and compare it with 
the expected entanglement of states chosen at random 
with respect to the Haar measure. We demand that the 
difference between the two expected entanglement values 
is consistently less than some small e. This will occur 
after some time t £ , for a given N and Na- Wc then check 
how t e scales with increasing TV with fixed Na, Fig- El 

These results suggest that the time required to achieve 
the Haar average is linear in the number of qubits and 
so we expect good scalability of our scheme [39f. It is 
worth noting here the relevance of the parameter <f>. Nu- 
merical indications suggest this is related to the abso- 
lute rate of convergence, with a maximum obtained by a 
value of approximately <j) = 5ir/8. It can be easily seen 
that two cases for which our scheme does not work are 
cf> = it and 4> = 2ir. These lead, respectively, to 2D and 
ID cluster state resources, for which fixed measurements 
are not universal. Interestingly, however, we have seen 
no evidence of any other general restrictions on permis- 
sible values of <fr. An intuitive explanation for this is 
that the set of possible single time-step evolutions, eq. 
([5]). are non-commuting operators. Taken together these 
may explore the space of unitaries with some complicated 
dependence on the parameter <$>. 
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FIG. 4: (Color online) Comparisons of the entanglement dis- 
tributions between the weighted graph states and states cho- 
sen at random with respect to the Haar measure, for N — 8 
and Na — 1,2,3,4 (plots from top to bottom). In both 
cases we calculate 10 5 samples and group into 500 evenly 
spaced bins. The distributions are then re-scaled such that 
J p(Sa) = 1. Black dashed lines represent the weighted graph 
state entanglements, and blue solid lines represent the Haar 
distribution. For comparison, we also include the distribu- 
tions for stabilizer states chosen at random [5] (solid lines 
capped with circles). 



C. Realization 

There are several methods and technologies in which 
states such as the one presented in Fig. [5]may be realized. 
Here we briefly mention two. Firstly, in optical lattices 
it is possible to implement two-body collisions that im- 



plement a Hamiltonian of the form H = Jo^P ® . To 
obtain the state in Fig. [2] it is then necessary and possi- 
ble to first implement this Hamiltonian between horizon- 
tal neighbors for a time t = it /(A J) followed by the im- 
plementation of the same pairwise Hamiltonian between 
vertical pairs for a time t = 4>/(2J) [3(|. The ability 
to generate the desired quantum state in parallel is con- 
trasted by the slight disadvantage that local addressabil- 
ity in such systems is difficult. Recently proposed s pin 
and polariton systems in arrays of optical cavities [3 111 
allows, in principle, for mechanisms analogous to the one 
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FIG. 5: (Color online) The scaling with N of the time taken 
for the average entanglement to agree with that associated 
with the Haar measure to a fixed accuracy e < 0.0001 (red 
squares), s < 0.001 (blue circles) and e < 0.01 (green trian- 
gles). For N = 2,3, 4 averages are over 10 6 realisations. For 
N = 5, 6 averages are over 10 5 realisations and for N = 7,8 
averages are over 10 4 realisations. In all cases 4> — 5ty/8. 

described above [32j but where individual nodes of the 
cluster state reside in distinct cavities and are therefore 
addressable. 

A somewhat different approach to achieve states such 
as those in Fig. [5] is via measurement induced interac- 
tion between cavities. The basic idea relies on the insight 
that one may first entangle the electronic degree of free- 
dom of atoms inside two distinct cavities with the cavity 
photons. These photons will then leak out of the cavi- 
ties where they are mixed on a 50/50 beamsplitter and 
then detected at the two outputs of the beam-splitter. 
This detection may implement a Bell projection on the 
atoms enabling the generation of entangled states [33| . 
Such Bell projections can be used to build large cluster 
states. These Bell projections may then actually be used 
directly to implement Controlled NOT gates [34[ in a loss 
tolerant way. Weighted graph states can then be gener- 
ated via the application of controlled phase gates with 
rotation angle <\>. These gates may be obtained from two 
controlled NOT gates supplemented by local rotations 



state by comparing the entanglement distribution of the 
output to that associated with the uniform distribution 
on states. Numerical results strongly suggest these two 
are identical, and so we may be conclude that the scheme 
is indeed very effective. 

From the theoretical perspective several important and 
interesting questions appear. Firstly one should prove 
rigorously that the scheme can generate any unitary op- 
eration. Furthermore, it may also be viewed be an inter- 
esting setting to study thcrmalization, since there is no 
classical randomness inserted by hand, yet the evolution 
appears to maximise the entropy of the input state. 
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V. CONCLUSION 



We have proposed a scheme to generate random quan- 
tum circuits. In our method, an experimenter prepares a 
particular type of very entangled state called a weighted 
graph state. The qubits carrying the input state are 
then entangled with part of the weighted graph state and 
the experimenter then just performs projective measure- 
ments in the |+)/|— ) basis. The randomness of the mea- 
surement outcomes chooses the particular circuit and no 
classical random numbers are necessary. We have tested 
the effectiveness of this scheme in randomising the input 
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APPENDIX A: GENERALISING THE SCHEME 

One could also consider generalising the scheme pre- 
sented in section IIII1 by including resource qubits in 
the state \ip) rjh = 7jfc|0) + 6 jk \l), along with gen- 
eral two qubit operations between columns. Further- 



more, we could also allow measurements in the ba- 
sis {(|0) ± e t6jk \l))/^/2)} . Defining new operators 

U [ l = (0\ a U ab \0) a + (-l) s e-^(l| a £/ ab |0) Q and lf$ = 
(0\ a U ab \l) a + (-l) s e- i; ^(l| Q t/ afc |l) a , a = r jk , b = r jk+1 
we can write the more general evolution operators as 



tff ?(ih_ h(0\bU L ] \0) b +S(0\ b U L ' \l) b 7 <0|b^ 10)6+5(016^ |1)6 
lVlj(0 ) — _q( k 1 _q<- k ~> _<?e=> _s< fc > 

\7<W \0) b + 6(l\ b U L > |l) b 7UW |0)6+<5(l|6C/ fl 3 \l) b J 



For operators U with entries u%j this can be written as, 
MAS k )= ( mn mi2 V (Al) 

Jy \ 77121 "722 J 

where 

mn = l(un + (-l) 5 ^ e~ l6]k u S i) 

+ 5{ui2 + {-l) S ^ e- i9ik u 3 2), 
mi2 = 7(iti 3 + (~l) Sj e- l9 > k u 33 ) 

+ 6(u u + (-l) s ^e- i0 ^ U34 ), 
m 2 i = 7(1*2! + (-l) S J* } e-^*U4i) 

+ S(u 22 + (-l) S ^e-^u 4 2), 
m 22 = 7(^23 + (-l^V^^s) 

+ 6{u 24 + (-l) s i k) e- lf> i k u A4 ). 



If t/ is a diagonal unitary operator, then we have the 
following restriction for unitary evolution, 

| 7 | 2 <iM33 + \S\ 2 u* 22 u A4 = 0. (A2) 

For some applications, such as hiding information about 
an input state, we are interested in obtaining a randomis- 
ing operation itself, as opposed to just a state with an 
expected amount of entanglement. This operation should 
ideally be unitary, and so the simplest guarantee is to de- 
mand that each step in the circuit also corresponds to a 
unitary operation. We then must impose the above con- 
dition from the outset, and so if we wish to adhere to the 
weighted graph state formalism {!]), then we are forced to 
apply only CZ gates between the neighbouring columns 

of Fig. m 



[1] D. P. DiVincenzo, D. W. Leung and B. M. Terhal, Quan- 
tum data hiding, IEEE Trans. Inf. Theory, 48(3) 580 
(2002). 

[2] J. Emerson, R. Alicki and K. Zyczkowski, Scalable Noise 
Estimation with Random Unitary Operators, J. Opt. B: 
Quantum Semiclass. Opt. 7, S347 (2005). 

[3] A. Scrafini, O.C.O. Dahlsten and M.B. Plenio, Thermo- 
dynamical state space measure and typical entanglement 
of pure Gaussian states, Phys. Rev. Lett. 98, 170501 
(2007). 

[4] A. Serafmi, O.C.O. Dahlsten, D. Gross and M.B. Plenio, 

Canonical and micro-canonical typical entanglement of 

continuous variable systems , J. Phys. A: Math. Theor. 

40, 9551 (2007). 
[5] J. Calsamiglia, L. Hartmann, W. Dur and H.-J. Briegel, 

Entanglement and decoherence in spin gases, Phys. Rev. 

Lett. 95, 180502 (2005). 
[6] J. Gemmer, A. Otte and G. Mahler, Quantum Approach 

to a Derivation of the Second Law of Thermodynamics, 

Phys. Rev. Lett. 86, 1927 (2001). 
[7] S. Popescu, A. Short, and A. Winter, The foundations 



of statistical mechanics from entanglement: Individual 
states vs. averages, Nature Physics 2, 754 (2006). 
[8] R. Oliveira, O. C. O. Dahlsten and M. B. Plenio, Efficient 
Generation of Generic Entanglement, Phys. Rev. Lett. 
98, 130502 (2007). 
[9] O. C. O. Dahlsten, R. Oliveira and M. B. Plenio, Emer- 
gence of typical entanglement in two-party random pro- 
cesses, J. Phys. A: Math. Theor. 40, 8081, (2007). 

[10] P. Hayden and J. Preskill, Black holes as mirrors: quan- 
tum information in random subsystems, Journal of High 
Energy Physics 09 (120) (2007). 

[11] A. Harrow, P. Hayden and D. Leung, Superdense coding 
of quantum states, Phys. Rev. Lett. 92, 187901 (2004). 

[12] J. Emerson, E. Livine and S. Lloyd, Random Circuits 
and Pseudo-Random Unitary operators for Quantum In- 
formation Processing, Science 302, 2098, (2003). 

[13] W. G. Brown, Y. S. Weinstein and L. Viola, Quantum 
pseudo-randomness from cluster-state quantum compu- 
tation, Phys. Rev. A 77, 040303(R) (2008). 

[14] J. Emerson, E. Livine and S. Lloyd, Convergence con- 
ditions for random quantum circuits, Phys. Rev. A 72, 



8 



060302(R) (2005). 

O. C. O. Dahlsten, Typical Entanglement-from the ab- 
stract to the physical, PhD Thesis, Imperial College [32] 
(2008). 

P. Hayden, D. W. Leung and A. Winter, Aspects of 
generic entanglement, Comm. Math. Phys. 265(1), 95 [33] 
(2006). 

E. Lubkin, Entropy of an n-system from its correlation 
with a k-reservoir, J. Math. Phys 19(5), 1028 (1978). [34] 
S. Lloyd and H. Pagels, Complexity as Thermodynamic 
Depth., Ann. of Phys. 188, 186 (1988). 
D. N. Page, Average entropy of a subsystem, Phys. Rev. 
Lett. 71, 1291 (1993). [35] 
S. K. Foong and S. Kanno, Proof of Pages conjecture on 
the average entropy of a subsystem, Phys. Rev. Lett. 72, 
1148 (1994). [36] 
A. Harrow and R. Lowe, Random Quantum Circuits are [37] 
Approximate 2-designs, quant-ph/0802.1919vl. [38] 
M. Hein, W. Diir, J. Eisert, R. Raussendorf, M. Van den 
Nest and H. J. Briegel, Entanglement in Graph States 
and its Applications, E-print arXiv:quant-ph/0602096vl; 
D. E. Browne and H. J. Briegel, One-way Quan- 
tum Computation - a tutorial introduction, E-print 
|arXiv:quant-ph/06032"26y i. 

S. Anders, M. B. Plenio, W. Diir, F. Verstraete and H. J. 
Briegel, Ground state approximation for strongly inter- 
acting systems in arbitrary dimension, Phys. Rev. Lett. 
97, 107206 (2006). 

L. Hartmann, J. Calsamiglia, W. Diir and H. J. Briegel, 
Weighted graph states and applications to spin chains, 
lattices and gases, J. Phys. B: At. Mol. Opt. Phys. 40 
S1-S44 (2007). 

D. Gross and J. Eisert, Novel Schemes for Measurement- 
Based Quantum Computation, Phys. Rev. Lett. 98, 
220503 (2007). 

D. A. Meyer and N. R. Wallach, Global entanglement in 
multiparticle systems, J. Math. Phys. 43, 4273 (2002). 
G. Smith, D. W. Leung, Typical Entanglement of Stabi- 
lizer States, Phys. Rev. A 74, 062314 (2006). 
O. C. O. Dahlsten and M. B. Plenio, Entanglement prob- 
ability distribution of bi-partite randomised stabilizer [39] 
states, Quant. Inf. Comp. 6, 527 (2006). 
D. Deutsch, A. Barenco and A. Ekert, Universality in 
Quantum Computation, Proc. R. Soc. Lond. A 449, 669 
(1995). 

O. Mandel, M. Greiner, A. Widera, T. Rom, T. Hansen 
and I. Bloch, Controlled collisions for multi-particle en- 
tanglement of optically trapped atoms, Nature 425, 937 
(2003). 

M. J. Hartmann, F. G. S. L. Brandao and M. B. Plenio, 
Nature Phys. 2, 849 (2006); A. D. Greentree, C. Tahan, 
J. H. Cole and L. C. L. Hollenberg, Nature Phys. 2, 856 



(2006); D. G. Angelakis, M. F. Santos and S. Bose, Phys. 
Rev. A 76, 031805(R) (2007). 

M. J. Hartmann, F. G. S. L. Brandao and M. B. Plenio, 
Phys. Rev. Lett. 99, 160501 (2007); D. G. Angelakis and 
A. Kay, New J. Phys. 10, 023012 (2008). 
S. Bose, P. L. Knight, M. B. Plenio and V. Vedral, Phys. 
Rev. Lett. 83, 5158 (1999); D. E. Browne, M. B. Plenio 
and S. F. Huelga, Phys. Rev. Lett. 91, 067901 (2003). 
I. E. Protsenko, G. Reymond, N. Schlosser, and P. Grang- 
ier, Phys. Rev. A 66, 062306 (2002); Y. L. Lim, S. D. 
Barrett, A. Beige, P. Kok and L. C. Kwek, Phys. Rev. A 
73, 012304 (2006). 

A. Barenco, C.H. Bennett, R. Cleve, D.P. DiVincenzo, 
N. Margolus, P. Shor, T. Sleator, J. Smolin and H. We- 
infurter, Phys. Rev. A 52, 3457 (1995). 
ff:|0/l)^| + /-). 

CZ: 1 00) h-> |00), |01) i > |01), 1 10) i — * |10), |11) i-> -}11). 
One could attempt a proof along the lines of [23 ]. by 
showing that the operators effected by the various mea- 
surement outcomes can be strung together to generate 
any linear combination of elements spanning the Lie Al- 
gebra u(2 N ). However, we have now restricted ourselves 
to specific one-parameter families of unitary operators, 
and so can no longer invoke the same generic arguments 
conventionally applied in such approaches. This leads to 
many difficulties. For example, when trying to demon- 
strate the appropriate irrationality conditions for the 
eigenvalues of the evolution operators explicitly, the re- 
sulting analytic expressions quickly become intractable. 
Thus, we have so far been unable to prove that fixed mea- 
surements and weighted graph state resources generate a 
universal set. It should be mentioned, however, that one 
can immediately see why fixed measurements on graph 
state resources are not universal. This would correspond 
to a parameter value of <j> = 7r in our model, for which the 
eigenvalues of the resulting unitary operators each con- 
tain one of the pairs ±7r/3 or ±27r/3. As these operators 
are periodic, we can not generate arbitrary linear com- 
binations of their algebra elements, and so we can never 
cover the full space of unitary operations. 
Note that the 'time' here refers to the number of columns 
in the weighted graph state that have been measured. 
Care should therefore be taken in comparing it with the 
number of two qubit gates necessary in a circuit model 
randomization. Note that for each time step here, 2N — 1 
two-qubit gates are needed, which effectively performs 
N — 1 two-qubit gates on the input state, even though 
only local measurements are made. The experimenter 
however needs to perform successfully (2iV — 1)Z two qubit 
gates in total, where I is the length of the circuit to be 
implemented. 



